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ABSTRACT 

We have determined the metallicity (0/H) and nitrogen abundance (N/0) of 
a sample of 122751 Star Forming Galaxies (SFGs) from the Data Release 7 of the 
Sloan Digital Sky Survey (SDSS). For all these galaxies we have also determined 
their morphology and obtained a comprehensive picture of their Star Formation 
History (SFH) using the spectral synthesis code STARLIGHT. The comparison 
of the chemical abundance with the SFH allows us to describe the chemical 
evolution of the SFGs in the nearby universe (z < 0.25) in a manner which is 
consistent with the formation of their stellar populations and morphologies. 

A high fraction (45%) of the SFGs in our sample show an excess of abundance 
in nitrogen relative to their metallicity. We also find this excess to be accompa- 
nied by a deficiency of oxygen, which suggests that this could be the result of 
effective starburst winds. However, we find no difference in the mode of star for- 
mation of the nitrogen rich and nitrogen poor SFGs. Our analysis suggests they 
all form their stars through a succession of bursts of star formation extended over 
a few Gyr period. What produces the chemical differences between these galaxies 
seems therefore to be the intensity of the bursts: the galaxies with an excess of 
nitrogen are those that are presently experiencing more intense bursts, or have 
experienced more intense bursts in their past. We also find evidence relating the 
chemical evolution process to the formation of the galaxies: the galaxies with 
an excess of nitrogen are more massive, have more massive bulges and earlier 
morphologies than those showing no excess. Contrary to expectation, we find 
no evidence that the starburst wind efficiency decreases with the mass of the 
galaxies. As a possible explanation we propose that the lost of metals consis- 
tent with starburst winds took place during the formation of the galaxies, when 
their potential wells were still building up, and consequently were weaker than 
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today, making starburst winds more efficient and independent of the final mass 
of the galaxies. In good agreement with this interpretation, we also find evidence 
consistent with downsizing, according to which the more massive SFGs formed 
before the less massive ones. 

Subject headings: Galaxies: stellar content — Galaxies: abundances — Galaxies: 
starburst — Galaxies: evolution — Galaxies: formation 



-4- 



Introduction 



One of the most important achievements of modern astronomy is the discovery of 
the process of nucleo-synthesis of chemical elements in stars. Once properly understood 
this process, coupled with the concept of stellar evolution regulated by the different 
masses of stars, gives us a unique insight about the chemical evolution of galaxi es 



(Edmunds fc Page 



1978 



Brodie fc Huchra 



1991 



Zaritskv et al. 



1994; 



Coziol et al. 



1998 



1999 



Henry et al. 



2000 



Pilyugin et al. 



2003 



Pilvugin & Thuanll201ll : 



Torres-Papaqui et al 



201 ll ). For instance, it is now well accepted that oxygen and sulfur are two elements 



produ ced by massive stars ( M > while 



stars ( IRenzini &: Voli 



1981 



McCall et al 



n itrogen is mostly a product of 



1985 



Evans fc Dopita 



1985 



Garnett 



ower mass 



1990)). 



Consequently, due to the longer time passed on the main sequence by stars with decreasing 



(Matteucci & Tosi 


1985; 


Garnett 


1990; 


Molla et al. 



2006 



Richer fc McCall 



20081 ). Assuming the initial mass function (IMF) does not vary 
between spiral galaxies, such time delay, when properly documented, may thus reveal 
something fundamental about how these systems formed their stars; for example, allowing 
to disti nguish between constant star formation over the formation by a se quence of stellar 



bursts ( Lehnert fc Heckman 



1996 



Coziol et al. 



1999 



Tremonti et al. 



20041 ) . 



To determine the "normal" or standard chemical evolution of Star Forming Galaxies 
(SFGs) in the nearby universe (z < 0.25), the observation and study of a large, homogeneous 
and statistically representative sample is req uired. This is where the Sloan Digital Sky 



Surve y (SDSS) project becomes so valuable ( lYork et al. 



2000 



Hogg et al. 



2001 



Pier et al. 



20031 ). By applying different automatic algorithms to the enormous data bank produced 
by the SDSS it is now possible to retrieve one of the largest and homogenous sets of 
spectral line ratios necessary to estimate the basic chemical abundances, and to describe 
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the chemical evolution process of SFGs in a manner consistent with their star formation 
histories and morphologies. 



Yin. Liang fc Zhang 



In recent articles (ITremonti et al.l 12004 ; iNagao et al.l l2006t llzotov et al.l 12006 



20071 ) data from the SDSS were already used to verify the consistency 



of the different methods devised in the past to determine the abundance of elements. One 
of the difficulties encountered in these studies is related to the rarity of the [OIII]A4363 
line. Theoretically, this line was recognized as crucial in order to determine an accurate 
temperature for the gas in HII regions. Unfortunately, [OIII]A4363 can only be observed 
in very low metallicity SFGs, which form a minute fraction of the SDSS galaxies. This 
fact emphasizes the importance of developing different empirical method s like R?x or R% to 



obtain chemical information 



McCall 



1984; 



Vacca fc Conti 



or a significantly larger sample of galaxies ( 



1992 



Thurston et al. 



1996 



Thuan et al. 



Pagel et al. 



2010) 



19791 : 



Recently, some authors (e.g. 



Nagao et al. 



2006 



Yin. Liang &: Zhang 



20071 ) advocated 



that we need to modify some of the empirical methods applied in previous abundance 
studies using only galaxies where [OIII]A4363 was observed. These authors based their 
claims on the fact that they found apparent significant differences in the abundances 
determined when they use their new calibrations. However, such point of view is somewhat 
problematic, as it assumes all galaxies follow the same chemical evolutionary pattern, 
independently of their mass or morphology, and assumes no evolution with redshift, while 
these are two ass umptions that need to be verified separa tely. Moreover, other researchers 



in the field (e.g. 



Tremonti et al. 



20oi 



Izotov et al. 



20061 ) that also tested thoroughly the 



empirical methods using SDSS data, have demonstrated that in general new empirical 
relations show results that are in good agreement with wha t was found before. For these 



reasons, but also for comparison sake with what was done in lCoziol et al. 



(119991 ). we choose 



for our study to apply the same empirical relations that were used before, but limiting our 
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chemical study to the t wo m ost important abundance ratios, O/H and N/O, which were 



shown by 



Izotov et al. 



(120061 ) to be less dependent on the method adopted to determined 



the gas temperature. 

Another important difficulty encountered in chemical evolution studies of SFGs is 
related with the contamination of emission lines by absorption features produced by 
the underlying older stellar populations. In our research we have solved this problem 
by subtracting a stellar population te mplate from each spectrum , as determined by the 



spectral synthesis code STARLIGHT (ICid Fernandes et al. 



20051 ). This method also has 



the advantage that through the fitted templates the star formation histories (SFHs) of the 
SFGs can be deduced, and other metallicity-independent parameters like the stellar velocity 
dispersion, that combined with the effective radius can yield an estimate of the mass of the 
bulge. In our study we use these new information in parallel with the morphologies which 
were determined independently to complete our view about the chemical evolution and 
formation process of SFGs in the nearby universe. 

This study is organized in the following manner. In Section [2] we describe how our 
sample of SFGs was constructed and how the data for our analysis were obtained. In 
Section [3] we present our results for the chemical abundances, and show how they varied 
with the mass and morphology of the galaxies. In the same section, we also explore the 
relation between the chemical abundances and the SFHs, and compare our results using 
STARLIGHT with some relevant models from Starburst 99. In Section H] we discuss 
our observations and propose a new interpretation. Our main conclusion can be find in 
Section |5j Many of our results were verified using statistical tests, which were regrouped in 
Appendix A. 
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2. Determination of the samples and data for analysis 



The data for our study were taken from the SDSS Data Release 7 (lAbazajian et al 



20091 ). Using the STARLIGHT Virtual Observatory servicqj we retrieved the spectroscopic 
data for 122751 SFGs, with emission lines having a signal to noise ratio S/N > 3 (adjacent 
continuum S/N > 10) and a redshift z < 0.25. The spectra were correc ted for Galactic dust 
extin ction and processed through STARLIGHT spectral synthesis code ( jCid Fernandes et al. 
20051 ). producing for each galaxy a stellar population template-corrected spectrum from 



which emission line fluxes were measured automaticall y. Note that STA RLIGHT estimates 



20071 ). which means 



and applies a correction for internal dust extinction ( lAsari et al. 
that the stellar population templates, and the star formation history deduced from these 
templates, as well as the template-subtracted emission spectra are free of this effect. 

In SDSS, the fibers have a fixed aperture and are centered on the nucleus of the 
galaxies. Comparing the physical projection of this aperture on the sky with the efficient 



radius as determined by 



Simard et al 



( 120111 ). we found that, except at very low redshifts, 
the spectra always cover the same physical regions comparable with the sizes of the bulges 
of the galaxies. This implies that our stellar population study and chemical analysis are 



mostly concerned with processes affecting the bulges of the galaxies. 



2.1. Separation between SFGs and AGNs 



The SFG c 



( iBaldwin. et al. 



assification presented in Figure [TJ is based on one diagnostic diagram 



1981 



Veilleux fc Osterbrock 



1987) , where we hay e 



criterion between AGNs and SFGs proposed by iKauffmann et al. 



standard interpretation of stellar ionized HII regions (iMcCall et al 



app 



(2003). 



ied the separation 
According to the 



1985 



Evans fc Dopita 



1 http:/ /www. starlight. ufsc.br 
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19851 ; ICoziollll996l ). the SFGs trace a sequence of increasing metallicity as the emission ratio 
[OIII]A5007/H/3 decreases, which is produced by the cooling effect of oxygen. Similarly, 
the abundance of nitroge n is also expected to grow in the same d irection, as the ratio 



[NII]A6584/Ha increases (IThurston et al. 



1996 



van Zee et al. 



1998|). 



Coziol et al. 



(119991 ) have determined 



Consistent with the standard interpretation, 
that some SFGs, especially those experiencing a starburst in their nuclear regions, present 
an excess of nitrogen abundance compared to the expected chemical evolution sequence 
traced by HII regions in late-type spirals. This is ill ustrated in Figure BK , where we show 



Coziol et al 



(119991 ): for the same 



the results of two ionization models as produced by 
metallicity, galaxies following the Sequence 2 (Seq. 2) have a nitrogen abundance about 0.2 
dex higher than galaxies following the Sequence 1 (Seq. 1). Figure [2] reveals that most of 
the SFGs in our sample have a data point that falls between these two models. Using the 
median of these two models, we have thus separated our SFG sample in two groups (see 
Figured). The SFGs that are nitrogen poor are identified as SFG 1, and represent 55.2% 
of the sample, while those that are nitrogen rich are identified as SFG 2, and represent 
44.8%. 

In Figure [2^ we notice that the path of Seq. 2 seem s to form an upper limit rather 



than a median as was determined by 



Coziol et al 



(1l999l ). The reason for this difference 



is because in 1999 the concept of transition objects (TOs), which are defined roughly as 
galaxies where both an AGN and intense star forming activity are present, was not applied 
systematically in classification studies, and the limit between SFGs and AGNs w as then 
located at higher values of the ratio [NII]A6584/Ha. Although 



Coziol et al. 



(119991 ) in their 



analysis presented evidence against the presence of an AGN in their sample, this possibility 
remains in this much larger and differently selected sample. Therefore, the SFG 2 galaxies 
may show an excess of excitation due to an AGN instead of an excess of nitrogen abundance. 



-9 - 



Interestingly, in Figure |2b we can see that the median of the two abundance models defined 
in (a) almost coincides with a different separation criterion between SFGs and AGNs which 



was previously proposed by 



Stasinska et al. 



( 120061 ). If we use this distinction we could 



separate our sample into two different new groups, according to which 64.8% of the galaxies 
are pure SFGs (SFG 3), and the rest (35.2%) are AGNs (SFG 4). 

For our sample of SFGs we have therefore identified two different interpretations for 
the increase of nitrogen emission in the standard diagnostic diagram. The first one suggests 
some SFGs show an excess in nitrogen abundance, the other suggests an excess in emission 
due to a supplementary source of ionizing photons, consistent with a central AGN. In order 
to discriminate which of these two hypothesis is the mos t probable, w e here present the 



Coziol et al. 



( 119991 ). consists in 



results of two tests. The first test, which was devised by 
comparing the ratio [NII]A6584/Ha with the ratio [SII]AA6717,6731/Ha:. The principle of 
this test is that if there is an excess of excitation, this excess should appear in both line 
ratios at the same time. 

We show the results of applying the first test to our SFG samples in Figure [3J The 
effect for an excess of excitation is not observed in the SFG 1 (a), nor than in the SFG 3 (c). 
Instead, we observe that at the same time the emission of nitrogen increases, the sulfur 
emission decreases. The same trend is also observed in Figure EJ^b) for the SFG 2, and in 
(d) for the SFG 4, although now most of the galaxies are well into the decrement phase of 
sulfur. Only an insignificant number of galaxies in these two samples are consistent with an 
excess of excitation in both lines. From this first test we conclude that, in general, we see 
no evidence for an excess in excitation that could have been produced by an AGN. 

As a second test, we present in Figure for the SFG 1 and SFG 2 and in Figure Hb 
for the SFG 3 and SFG 4 the graph of the Ha luminosity as a function of the luminosity of 
the continuum at 4800 A. The continuum fluxes were measured in the raw spectra, that is, 
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before they were treated by STARLIGHT, and were only c orrected for Galac tic extinction. 
The principle of this test is based on a suggestion made by lOsterbrockl (119891 ). according to 



who the linear relation (in Log) between the two luminosities must be different for galaxies 
excited by AGNs and galaxies excited by stars, because of the different powers laws in 
these objects relating the continuum with the ionizing flux. In Figure H] we observe that 
the linear relations for our four SFG groups have almost the same slopes. These linear 
relations have extremely high factors of correlation (see Tabled]). For comparison we have 
traced the relation found fo r two sample of SDSS AGNs classified as Seyfert 2 and LINERs 



( ITorres-Papaqui et al. 



20121 ). The slopes of the linear relations for the SFGs are significantly 
shallower than for those measured in AGNs. From this second test we conclude, once again, 
that there is no evidence for the presence an AGN in our sample of SFGs. 

Taking into account the high level of consistency of the above two tests, we conclude 
that in general we see no evidence for AGNs in our sample of SFGs, and we find 
consequently no reason to change o ur classification based on the separation criterion 



proposed by iKauffmann et al. 



( 120031 ) . For the remaining of our analysis we will keep only 



the distinction between nitrogen poor (SFG 1) and nitrogen rich (SFG 2) SFGs. 



2.2. Determination of abundances and physical parameters 



Once the possibility of an AGN as a source of ionization for the gas is eliminated, 
standard methods for HII regions can be applied to determine the chemi cal abundances 



of the SFGs. For the metallicity, O/H, we used the method described in 



Vacca fc Conti 



199 2J), which yields va 



( lEdmunds fc Pa ge! 



described in 



ues with a typical uncertainty of the order of ±0.2 dex 



1984). For th e nitrogen abundance, N/O, we utilized the method 



Thurston et al. 



(119961 ) . where the temperature of the gas, T/j, is first estimated 



using the ratio ( [Oil] AA3727, 3729 + [OIII]AA4959, 5007) /H/3, which is then combined with 
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the ratio [NII ] AA654 8, 6584/[OII1AA3727, 3729 to estimate the ratio N/O. According to 



McCall et al. 



(|1985[ ). the dispersions in the distributions of these two line ratios are due to 
a variation in metallicity (O/H) and relative nitrogen abundance (N/O), which justify the 
abundance determination method. A cor rection for ionizati on may become important only 



at low metallicity, 12 + log(0/H) < 8.2 (jlzotov et al 



20061), which is not the case for the 



galaxies in our sample. When necessary, all our abundanc es and results f r om m odels are 



given relative to the new Sun metallicity, as determined by 



Asplund et al. 



( 120041 ) 



From the STARLIGHT template we retrieved the mean age of the stellar populations, 
from which we can deduce the SFH of the galaxy. The SFH describes how the star 
formation rate (SFR) in a galaxy has varied over its lifetime. Since the stellar templates 
are corrected for internal extinction the SFR are also free of this effect. To produce the 
SFR we worked with a smoothed v ersion of the population vector obtained from the output 



of STARLIGHT flAsari et al. 



20071 ). where at each time, t iy the SFR(tj) is co mposed of a 



mixtu re of stars taken from six different metallicity groups. As discussed by 



Asari et al. 



(120071 ) in their section 5, any value in the range 25 Myr bin for the stellar library yields a 
strong correlation with the luminosity in Ha, which is used to measure the current SFR. 
However, this bin in age is not critical in STARLIGHT synthesis code, because values in the 
whole 10-100 Myr yield correlations of similar strength (this excludes WR stars which have 
much shorter ages, but there is no trace of such stars in our spectra). This range of ages is 
fully consistent with what is used in other synthesis codes, like for example Starburst 99. 
We therefore expect our results using STARLIGHT to be fully comparable and consistent 
with Starburst 99 models (see Section 3.3). 

Also from the STARLIGHT templates, we have retrieved the stellar vel ocity dispersio n 



of th e stars, a*. After correcting for the instrumental resolution (following iGreene fc Ho 
20061 ). the <7* are used in combination with the effective radii determined by 



Simard et al 
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( 120 111 ) to estimate the masses of the bulges (applying the virial theorem). Note that within 
the limits in redshift of our study the spatial projection of the aperture of the SDSS fiber 
varies from 1 kpc to 7 kpc, but at the same time the effective radii of the galaxies also 
increase, which implies that the fraction of the galaxy covered by the fiber at any redshift is 
roughly the same. Because the projected apertures are comparable with the effective radii 
of the galaxies where the kinematic of the bulge is dominant, no correction for the effect of 
rotation was applied to cr*. 

As an independent parameter (that is, not depending on STARLIGHT), we have 
classified all the galaxies in our sample adopting a morphological index T, varying on a scale 
from -5 to 10 (—5 = E, 10 = Irr). The Hubble morphological types are presented in Table [2j 
togeth er with the correspondence adopted between the de Vaucouleur's morphological 



index (jde Vaucouleurs et al 



19911 1 and our own. We have determined the morphology 



using the correlations between photometri c colors, the inverse con centr ation index and 



the m orphological types, as was found by 



Shimasaku et al 



(1200 lh and 



Fukugita et al. 



(120071 ). The photometric colors are u — g, g — r, r — i, and i — z, which are defined in the 
photometric system of the SDS^|. The inverse conc entration corresponds to R 50 (r) / Rgo(r) , 



that compares the Petrosian radii (jPretosian 



19761 ) containing 50% and 90% of the total 



flux in the r band. A K-correcti o n was applied to the Petrosian magnitudes using the code 



developed by 



Blanton fc Roweisl (120071 ). 



A last independent parameter useful for our an alysis is the abso 



which was retrieved from the SDSS Data Release 7 (lAbazajian et al 



ute m agnitude in B, 



20091 ). Being a good 



indicator of star 



Gyr (jCoziol 



brmation activity over relatively long periods of time, from 0.5 up to 6 



19961 ). the luminosities in B can also serve as a proxy for the total mass of the 



galaxies. Considering the low limit in redshift of our study, a simple cosmology was applied, 



2 http:/ /casjobs. sdss.org 
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adopting an Hubble constant Hq = 75 km s 1 Mpc 



Results 



3.1. Variation of abundances, and relation with mass and morphology 

In Figure |5] we trace the abundance of nitrogen as a function of the metallicity, 



distinguishing between nitrogen poor (SFG 1) and nitrogen rich (SFG 2 
on this figure we show the two abundance sequence models obtained by 



SFGs. Overplotte d 



Coziol et al. 



f |l999h . 



Also shown are the relations for the differ ent closed-box chemical e v olutio nary models for 



the production of nitrogen as proposed by 



Vila- Cost as fc Edmunds! (119931 ): the secondary 



and the primary + secondary models. In general, the SFG 1 seem to follow the secondary 
enrichment model. This is obviously not the case for the SFG 2, however, which for the 
same value in log (O/H) show an excess of nitrogen, on average of 0.15 dex, compared to 
the SFG 1. 

In Figure |5b we also note a remarkable feature: the chemical enrichment process for the 
SFG 2 is not continuous-we observe an abrupt jump by 0.3 dex in nitrogen at a metallicity 
log (O/H) = —3.4 (almost the solar metallicity). The reason why we note this p articularity 



Coziol et al. 



Jl999j) for their 



is because exactly the same observation was done before by 
comparatively small sample of about 100 UV bright SBNGs. The fact that exactly the same 
feature appears in our significantly larger sample and generally selected sample (no other 
criterion other than the SFG classification was applied) is quite remarkable. It suggests 
that this feature must be a characteristic of SFGs (it affects 45% of all the galaxies in our 
sample), and not a peculiar trait of starburst galaxies. This phenomenon must be related 
consequently to some common mechanism or event occurring during the chemical evolution 
of all these galaxies. 
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In Figure [5^, although the nitrogen distribution for the SFG 1 fits the secondary 
relation, we may detect a slight trend toward relatively high nitrogen abundance also in 
this sample. This is better perceived in Figure \5jp where we have plotted over the data 
th e medians of the di fference in nitrogen abundance relative to the secondary relation. 



In 



Cozioletal 



( 119991 ) the UV bright starburst galaxies were showing a trend toward a 
decrease in the excess of nitrogen at high metallicity. The same trend may also be visible 
in Figure [5H for the SFG 2, although, fr om the medians, we conclude that the decrement is 



much weaker than what was observed in 



Coziol et al. 



(1999) 



In Figure M we trace the metallicity and nitrogen offset from the secondary relation 
as a function of the absolute magnitude in B. I n Figure and b, we have included the 



Tremonti et al. 



( 12004] ) for their sample of 



luminosity-metallicity relation as determined by 
53400 SDSS SFGs (correcting for the Hubble constant used in our study). Our samples 
taken as a whole (adding the SFG 1 with the SFG 2) are in good agreement with this 
relation. However, when we introduce the distinction between nitrogen rich and nitrogen 
poor SFGs we obtain a different perspective. The metallicity of the SFG 1 increases with 



Tremonti et all (120041 ). while it increases 



the mass more rapidly than the relation found by 
less rapidly in the SFG 2. In fact, and as the correlation tests in Table [3] are confirming, 
the metallicity of the SFG 2 stays almost constant. Therefore, it is like the SFG 2 become 
deficient in oxygen as the mass increases. This is verified in Figure [TJ which shows the 
different histograms for the metallicities and absolute magnitudes. Although the SFG 2 are 
on average more massive (more luminous in B) than the SFG 1, their peak in metallicity 
lays at a lower value than in the SFG 1. This confirms that the SFG 2 are relatively 
deficient in oxygen as compared to the SFG 1. Note that the same phenomenon was also 



observed by 



Coziol et al. 



(119971 ) for th e same samp 



were shown later to be nitrogen rich ( 1 Coziol et al. 



e of UV bright starburst galaxies that 



19991 ). So it seems now that the two 



phenomena must be somehow connected. 
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In Figure [6fc we find that the excess of nitrogen in the SFG 1 is slightly increasing 
with the mass (although in Table [3] we have a 36-38% percent chance that the correlation 
is spurious). On the other hand, in Figure IHli for the SFG 2 the excess stays constant 
independently of the mass. This is confirmed by the correlation tests (Table [3]). It seems 
therefore like the chemical abundances in the SFG 2 reach some sort of physical limit 
connected with the masses of the galaxies. 

In Figure [8] we compare the results obtained for the morphologies, the mean stellar 
population ages, the masses of the bulges, and absolute magnitudes in B, as measured 
in the SFG 1 and SFG 2 samples. The box-whisker diagrams show that the SFG 2 have 
earlier-types (Figure [8^) and older stellar populations than the SFG 1 (Figure |8b). The 
masses of the bulges in Figure [Ht are also increasing in the SFG 2, which is consistent with 
the earlier morphologies. The difference in absolute magnitude in B (Figure EH) also implies 
higher masses for the SFG 2 as compared to the SFG 1. All the above differences were 
confirmed at a level of confidence of 95% by statistical tests (see Figure [221 in appendix A). 
All these variations and differences point toward a connection between the chemical 
evolution process and the formation process of the SFGs. 



3.2. Analysis of SFH using STARLIGHT 

Through the fitting templates using STARLIGHT we have deduced the history of 
star formation (SFH) for all the galaxies in our sample. In order to investigate how 
the chemical evolution in the SFGs correlates with their SFHs, we have separated the 
diagram of abundances in Figure [9] in 4 subsamples with increasing metallicity. In each 
metallicity bin we further distinguish between low and high nitrogen abundance, obtaining 
in total 8 different subsamples. The corresponding SFH for these 8 subsamples are shown 
in Figure [101 The median physical characteristics of the galaxies in each subsample are 
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reported in Table 1 for the SFG 1 and Table [5] for the SFG 2. 

The SFH describes how the star formation rate (SFR) of a galaxy varied over its 
lifetime. In particular, the SFH allows to distinguish between galaxies that have experienced 
more or less constant star formation rates (CSFR), as is expected for late-type spirals, and 
galaxies in a present starburst phase (SBP). Galaxies with CSFRs have a flat SFH (e.g., 
subsamples 4 and 8 in Figure [TUl) while galaxies in SBP, with recent SFR higher than in its 
past, have a negative SFH slope (e.g., subsamples 1 and 5 in Figure [TO]) . Consequently, a 
galaxy with a positive SFH slope would have formed most of its stars in the past, as it is 
expected for early- type spirals and ellipticals (no case is observed in our sample). The SFH 
also yields information about the intensity of star formation: for any interval of time the 
SFR corresponds to the median, which implies the higher the median the more intense the 
star formation at this time; consequently, the higher the intensity of star formation and the 
higher the mass of stars formed, which is equal to the surface below the curve. 

For the SFG 1 sample in Figure UU] (left panels) we observe that almost all the galaxies 
with low nitrogen abundance (in subsamples 1, 2 and possibly 3) show some evidence of 
recent bursts of star formation. For the SFG 1 that are rich in nitrogen, we observe a 
transition to constant star formation starting with subsamples 7 (galaxies in subsamples 5 
and 6 are still in bursts). In the SFG 2 (right panels in Figure fTUll the sharp increase in 
nitrogen abundance seems connected with intense (more intense than in the SFG 1) bursts 
in star formation (subsamples 1, 5 and 6). We also note that the transition to constant 
star formation appears earlier than in the SFG 1, in subsample 2 instead of 3. In general, 
we observe that galaxies having higher intensity bursts or higher SFR at any time, end up 
with higher nitrogen abundance (comparing the continuous curves with the dashed ones). 
However, the difference seems to decrease at high metallicity (e.g., subsamples 4 and 8 in 
both the SFG 1 and SFG 2 samples). 
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The box-whisker plots for the SFG 1 in Figure [IT] and for the SFG 2 in Figure [12] 
are consistent with our previous analysis based on the abundances (see also Table [3J and 
Table H]). In both samples, we observe a change toward earlier-type morphology and older 
stellar populations as the metallicity and nitrogen increases. We also observe an increase 
in bulge mass and total mass. In general, the more massive and early-type SFGs with low 
metallicity are those that are experiencing more intense bursts, and at high metallicity 
the more massive and early-type SFGs are also those that have experienced higher level 
of star formation in the past. All these differences were confirmed by statistical tests in 
appendix A (see Figure 1231 and Figure UM . 



3.3. Comparison with Starburst 99 models 



The results of our analysis bas ed on STARLIG HT were compared using the model for 



star forming galaxies Starburst 99 (ILeitherer 



1999 



We used some of the original 1999 



dataset, considering only instantaneous or constant star formation scenarios, and have run 
a few more models using the WINDOWS version of the software. The main goal of this 
comparison is to verified if we can obtain Starburst 99 models that are consistent with our 
STARLIGHT analysis, in particular, by isolating the galaxies that show recent starbursts. 
Note that we did not performed an exhaustive study using Starburst 99, because this would 
have fall well beyond the limits of the present study. However, we have found that just a 
few simple modifications to the parameters of the original dataset were sufficient to yield 
results in good agreement with STARLIGHT. 

The only two parameters produced by Starburst 99 that can be compare with our 
results from STARLIGHT are the Ha equivalent width and absolute magnitude in B. The 



3 http:/ /www. stsci.edu/science/starburst99/docs/default. htm 
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first parameter is a good indicator for star formation over short timescales, for this is the 
ratio of the flux in emission with the flux in the continuum. As we already mentioned, 
the second parameter is a good indicator for star formation over longer timescales (0.5 to 
6 Gyr), and a reasonable proxy for the total mass of galaxies. In Figure [TBI we show the 
predictions of instantaneous bursts models, having different powers of the IMF (a) and 
upper mass limit (M up )] the Salpeter IMF corresponds to a = 2.35, and only two mass 
limits were tested, 30 and 100 M Q . The data presented are the medians and percentiles 
as measured in the galaxies found in the abundance subsamples defined in Figure [9] It 
is obvious that the absolute magnitudes of the SFGs are too high to be reproduced by 
an instantaneous burst scenario. This result clearly suggests that in these galaxies star 
formation happened over longer time scales. 

The next series of models in Figure [T4] are for constant star formation. The first 
good fits are obtained with a steeper IMF than Salpeter, a = 3.30, and upper mass limit 
M up = 100 M . But, then we loose most of the SFGs that STARLIGHT was identifying as 
starbursts. Based on this series of models it becomes obvious in which direction we need 
to modify the parameters to obtain fits that are consistent with STARLIGHT analysis. 
The IMF has to become steeper still, implying that there are less massive stars and more 
intermediate mass stars ionizing the gas in these galaxies (more B than O stars). But 
then to compensate, the total amount of ionizing photons produced has also to be higher, 
which justify the upper mass limit M up = 100 M . Our two best fits, which yield almost 
a one to one relation with the STARLIGHT results are presented in Figure [15j They have 
M up = 100 M Q and an IMF varying from a = 4.00 (a and d) to a = 4.10, (b and e). On the 
other hand none of these models reproduce the trend in metallicity we observed. 

Based our comparison, we do not claim we have determined that the IMF in the SFGs 
is much steeper than Salpeter or that the upper mass limit is 100 M . What we suggest, 
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rather, is that the results of Starburst 99 that are most consistent with STARLIGHT implies 
recent bursts of star formation happening over a background of constant star formation. 
But, considering that constant st ar formation in galaxies can be easily reproduced by a 



sequence of instantaneous bursts (ILeitherer 



19991 ). therefore, it seems that the best scenario 



that the models of Starburst 99, consistent with STARLIGHT, imposes on the data, can 
be, for both SFGs, a sequence of bursts of star formation over a re latively long pe r iod o 



time. Similar resu 
and 



Coziol et al 



ts were obtained before in starburst studies like 



Goldader et al 



(1997) 



(120011 ) . There is also now mounting evidence that typical starburst 



galaxies, like M 82, have experienced more than one burst spread both in ti me and space 



(Ide Griis et al. 



2006 



2001 



Mayya et al 



Smith fc Ga llagher 



2006 



2001 



Forster Schreiber et al. 



Strickland fc Heckman 



to starburst galaxies was already verified by 



2009). T 



2003a.b 



Smith et al. 



r at thi s scenario applies in general 



Coziol et al. 



(120011 ) . These last authors even 



confirmed the dominant presence of B type stars in the ionizing regions of th ese galaxies 



which is a natural consequence of a sequence of bursts over a few Gyr period (jCoziol 



199J). 



From our comparison with Starburst 99 models we conclude that our STARLIGHT 
analysis (except for the metallicity) yields results that are consistent in all the SFGs with a 
sequence of bursts of star formation over a long period of time (of Gyr scale). 



4. Discussion 



In 



Coziol et al. 



(119991 ) . the authors suggested that the discontinuous pattern observed 
in the diagram of nitrogen vs. metallicity was the result of a special mode of star formation 
in starburst galaxies. Assuming these galaxies experienced a succession of bursts over a 
few Gyr period, decreasing in intensity with time, and assuming a time delay between 
the production of oxyge n by massive O and B stars and the production of nitrogen by 



intermediate mass stars (IGarnett 



19901 ). then the higher the intensity of the first bursts, the 
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higher the increase of nitrogen compared to the oxygen. 



Although our present observations seem to be in good agreement with the above 
starburst scenario, our new analysis proposes a slightly diff erent interpretation . First, our 



Coziol et al. 



(119991 ) . The only 



sample of SFGs is more general than the sample defined in 
criterion applied in the selection of the galaxies is that they must not show any evidence 
of ionization from an AGN. Second, our analysis of the SFHs produced by STARLIGHT 
suggests that there is no dichotomy in the star forming mode of the SFGs, that is, all the 
SFGs follow the same star formation pattern, consistent with a sequence of bursts of star 
formation over a period of a few Gyr. The comparison of our STARLIGHT analysis with 
Starburst 99 models is extremely compelling on this point. What our new analysis suggests, 
therefore, is that the nitrogen excess is not an anomaly related to some special galaxies like 
starbursts, but is more likely a typical characteristic of SFGs in the nearby universe. This 
also implies that the phenomenon that creates the excess of nitrogen must be extremely 
common in star forming galaxies, although possibly more obvious in starburst galaxies. 

Consistent with the above interpretation, we believe that the common phenomenon 
responsible for the depletion of oxy gen and relative increase in nitrogen as observed in 



the SFGs could be starburst winds rtHeckman. Armus fc Miley 



1990 



Lehnert fc Heckman 



19961 ). In particular, this mechanism is suspected to be the main pheno menon responsib le 



for in creasing the abundance of metals in the intergalactic medium (see 



2005 



Veilleux et al 



, and refer ences therein ) . Now , it is easy to adapt the model of a sequence of bursts 



as proposed by 



Coziol et al 



(119991 ) including as a major component the effect of starburst 
winds. We assume that a starburst wind in a region is effective (that is, capable of ejecting 
metals out of the region) only for a short period of time, of the order of a few 10 7 yrs 
(as is confirmed using Starburst 99). Consequently the main effect of this wind would 
be to preferentially deplete the region of the chemical components produced by massive 
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stars (mostly oxygen, neon, iron, but also sulfur, argon and calcium; see Strickland & 
Heckman 2009, and references therein for a full discussion of starburst wind loads). Then, 
when intermediate stars begin producing nitrogen, a few 10 8 yrs after the beginning of 
the burst, the intensity of the wind would already be to low to have a significant impact 
on the abundance. Moreover, within the sequence of bursts model we do not expect the 
next generation of bursts to affect the nitrogen abundance in the region, simply because 
these bursts will happen in different regions, considering that in such a brief interval of 
time, a few 10 8 , the conditions in a region that suffered a burst would not be favorable for 
another one. However, these conditions could easily become favorable in another nearby 
region, consistent with the propagation of star formation hypothesis; after 10 8 yrs active 
star for mation regions would have propagated a few kpc from the first burst region (see for 



example 



Nomura fc Kamaya 



2001 



, and references therein). 



Can we verify the starburst wind hypothesis directly with our observations? To answ er 



Strickland fc Heckman! (J2009J), 



this question we refer to the detailed calculation made by 
according to who the difference between the nitrogen rich and nitrogen poor SFGs should 
depend on the intensity of the bursts. In their mod el, a starburst wind is cons i dered to b e 



1996 



Heckman 



effect ive above an intensity of 0.04 M Q yrs _1 kpc~ 2 (jLehnert fc Heckmanl 
20031 ). Using STARLIGHT we have estimated the intensity of the bursts (see Tabled]) 
for the galaxies in both samples showing the higher recent SFR (subsamples 1 and 5). In 
Figure [161 we show the corresponding box-whisker plots. We see that the intensity of the 
bursts in the galaxies with the most obvious excess of nitrogen (the SFG 2 in subsample 5) 
are significantly higher than in galaxies without an excess (S FG 1). Also shown in this figure 



is the threshold calculated by 



Strickland fc Heckmanl ( 120091 ). The median and average of 



the SFG 2 are well above this limit. Note that we should not expect a clear cut distinction, 
considering that these bursts are only the most recent of a series. However, statistically the 
difference is significative, and in good agreement with the starburst wind hypothesis. 



-22 - 



In Table [6] we also give the statistics for the absolute magnitudes in B. The SFG 2 
are more luminous in B than the SFG 1. This is consistent with other observations in 



starburst galaxies sugge sting that in general t 



the mass of the galaxy (IGoldader et al. 



1997 



re intensity of a s tarburst increases with 



Coziol et al. 



20011 ) . Also consistent with 



this difference of characteristic, the STARLIGHT analysis in Figure [TOl reveals that in the 
past the SFG 2 have also experienced more intense star forming episodes than the SFG 1. 
All these differences were found to be significant using statistical tests (see Figure [25] in 
appendix A). In general, our results are quantitatively consistent with the predictions made 
by the starburst wind model. 

On the other hand, one observation seems somewhat counterintuitive in regard to the 
starburst wind hypothesis. Many authors have discussed and concluded that the higher the 
mass of a galaxy and the more difficult it should be for this system to loose metals through 



starburst winds (e.g. 



Tremonti et al. 



20041 ). Yet, our observations show no evidence of such 
trend. In Figure [6H for the SFG 2 the excess of nitrogen stays constant independently 
of the mass, and in Figure |6b we may even detect an increase of nitrogen with the mass 
in the SFG 1. On this matter, our analysis of STARLIGHT results in parallel with our 
determination of the morphologies and mass of the bulges is quite revealing. It allows us 
to connect the excess of nitrogen with massive galaxies, having developed massive bulges 
and earlier type morphologies. At the same time it shows that the star formation in the 
past of these galaxies was more intense than in those without the nitrogen excess. There is 
consequently a very strong relation with the formation process of the galaxies. This may 
suggest that the abundance differences observed between the SFG 2 and SFG 1 took place 
mostly during the formation process of these galaxies. This could have happened when the 
gravitational potential wells of these galaxies were less massive, because still growing in 
size, making the intense starburst winds slightly more efficient and independent of the final, 
present mass. 
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Consistent with the formation hypothesis, it was recently reported by lPilyugin &: Thuan 



(120 111 ) that the SFGs are showing evidence for downsizing, which states that massive 
galaxies formed at higher redshifts than smaller mass galaxies. If we can find such evidence 
in our sample of galaxies, consistent with the scenario we proposed for the excess of 
nitrogen, then we would have established another link between the chemical evolution and 
the formation process of these galaxies. 

To test the downsizing hypothesis in our sample we have separated the SFG 1 and 
SFG 2 samples in 9 redshift subsamples increasing by 0.025. The different bin ranges are 
given in Table [7J for the SFG 1 and Table [H] for the SFG 2, together with the number 
of galaxies in each redsh ift subsamples. Note that our limit in redshift is lower than in 



Pilyugin fc Thuan! (120111 ). z= 0.25 instead of z= 0.4. In Figure [T71 we show the box-whisker 
plots for the variations of the properties of the SFG 1 as a function of the redshift. The same 
graph is shown for the SFG 2 in Figure [181 In the SFG 1 we see an increase in metallicity 
and nitrogen excess at higher redshift. This is not observed however in the SFG 2 sample, 
where the metallicity and nitrogen excess stay almost constant (see Figure [26] and Figure [27] 
in appendix A for confirmation by statistical tests). Note that this difference between the 
SFG 1 and SFG 2 goes against what would be expected if the variations observed were due 
to an aperture effect-the same effect should have been observed in both samples. Moreover, 
in both samples, we observe a comparable increase in bulge masses and absolute magnitudes 
(cf. Figure [26] and Figure [27] in appendix A). Because these two observations are clearly 
independent of the aperture, this result suggests the masses of the galaxies increase at 
higher re dshift. Taken at face valu e, therefore, our observations are comparable to those 



made by 



Pilyugin fc Thuan] (120111 ). suggesting that the most massive galaxies formed first 



at high redshift. 



However, we do not agree with the assumption made by 



Pilvugin fc Thuan! (1201 lh that 
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the difference in abundance implies a decrease of SFR at high redshift. For example, in 
Figure [T7J and Figure [IS] we observe an increase of late-type spirals and of younger stellar 
populations at higher redshift (cf. Figure |2"o1 and Figure [27] in appendix A). These trends 
seems somewhat in contradiction with the other observations described above. Indeed, how 
can the bulges increase in mass at high redshifts while the number of late-type spirals and 
number of young stellar population also increase? The most plausible explanation is that, 
in fact, the SFR in both samples increases at higher redshift. Confirming this interpretation 
we show in Figure [19] for the SFG 1 and Figure [20] for the SFG 2 the variation of SFH 
with redshift. For the SFG 1, all the galaxies at any redshift show a SFH consistent with 
constant SFR. There is consequently no evidence of a decrease of star formation with the 
redshift. Similarly, all the SFG 2 at any redshift show a starburst like nature (the peak at 
young age is higher than the peak at old age) . This difference is consistent with their excess 
of nitrogen: at any redshift the SFG 2 are experiencing more intense bursts than the SFG 1 
and consequently show more clearly the effect of starburst winds. However, and this is the 
most significant result yielded by STARLIGHT, in both samples we also observe a general 
increase of the SFR with the redshift. This is fully consistent with downsizing, showing 
that the more massive galaxies formed first at higher redshifts. 

Note that according to our analysis the downsizing phenomenon can be observed at 
almost any small increment of reshift. This was unexpected and somewhat surprising. Can 
we be wrong in our interpretation? That is, can this be, instead, the product of some sort 
of luminosity bias, the galaxies are getting bigger and more massive because they are more 
luminous at high redshift? In our sample we do observe a decrease in the number of SFGs 
with the redshift (cf. Tables 7 and Tables 8), which, taken at face value, may look as an 
argument in favor o f a luminosity bias. However, this is inconsistent with the study of 



Strauss et al. 



( 120021 ). where the criteria for the SDSS spectroscopic survey are presented 



and discussed, and the authors conclude: "The uniformity and completeness of the galaxy 
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sample make it ideal for studies of large-scale structure and the characteristics of the galaxy 
population in the local universe." (We emphasized the completeness.) 

In fact, the explanation for the "missing" galaxies is simple. Our sample is indeed 
biased, but in another way, which is that all the galaxies are SFGs, with a S/N in line ratios 
> 3. We verified that there are no missing galaxies in the SDSS survey at high redshift, 
simply they are not SFGs and do not satisfy our S/N criterion. In particular, we have 
found an important number of narrow line emission galaxies with weak or missing emission 
lines that complete the sample. All these galaxies have earlier morphological types than the 
SFGs and are found in clusters and groups (the majority are consistent with low luminosity 
AGNs). In other words, the SDSS survey shows that the activity type of galaxies changes 
at higher redshift, including less and less SFGs, which is consistent with downsizing. 

Within the hierarchical biased galaxy formation paradigm, assuming the starburst 
phenomenon is related with galaxy formation, the SFGs, which are all spirals, are typical of 
the field, which means they formed relatively recently from low density fluctuations. At this 
level of density fluctuations the differences observed between the SFG 1 and SFG 2 must 
be related to very local variations, that is, on small scales of mass density, those leading to 
the formation of massive bulges. The downsizing phenomenon emphasizes that the redshift 
is really cosmological in nature. Rather than a variation of distance, it really describes a 
variation of state of the universe, namely here its density. 



5. Conclusion 



In ICoziol et al.l (119971 1. it was reported that U V bright star b urst g alaxies are deficient 
in oxygen compared to normal late-type spirals. In ICoziol et al.l ( 119991 ) it was also reported 
that the same sample of starburst galaxies show a possible excess in nitrogen abundance. In 
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the present study, using a much larger (122751) and generally selected sample of SFGs (that 
is, the only restriction was that they do not show the presence of an AGN), we find exactly 
the same phenomena, suggesting that this is a common trait of star forming galaxies, and 
not a peculiar characteristics of starburst galaxies. 

Our analysis of the SFH using STARLIGHT, in good agreement with Starburst 99 



models, suggests that the depletion in oxygen and the relative excess in nitrogen are m ost 



probably due to the effect 


of intense star 


mrst winds ( 


Hec 


mian, Armus & Milev 


1990; 


Lehnert & Heckman 


1996: 


Heckman 


2003; 


Tremonti et al. 


2004 




Strickland & Heckman 



20091 ) happening during a prolonged sequence of bursts of star formation. According to our 



analysis, all the SFGs form their stars through a sequence of bursts and what produces the 
chemical differences is a variation in intensity of the bursts. The SFGs experiencing the 
more intense bursts (with a median intensity above a particular threshold), or those that 
have experienced more intense bursts in the past, all show the effect of starburst winds: 
they are deficient in oxygen and relatively rich in nitrogen (N/O). We illustrate our model 
in Figure I2T1 

Contrary to expectation, however, we find no evidence that the efficiency of the 
starburst wind decreases with the mass of the galaxy. Instead, our data suggests the 
intensity of the bursts grows with the mass, and with it the depletion of oxygen and 
the excess of nitrogen. To explain this observation, we propose that the abundance 
characteristics of these galaxies took shape during their formation process, when their 
gravitational wells were still forming, and consequently less massive and susceptible to 
decrease the efficiency of starburst winds. Consistent with this interpretation we have 
shown that the galaxies with an excess of nitrogen are the more massive, have bigger bulges 
and have an earlier morphological type than those without an excess. We have also found 
evidence consistent with downsizing, suggesting that the most massive galaxies formed first 
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at higher redshift. 

Considering the generality of our analysis we conclude that the formation process of 
the SFGs is an open process, the galaxies loosing mass and energies to their environment. 
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A. Results of statistical tests 



The statistical framework of the tests we used is based on a new parametric 



AN 
by 



OVA model introduce d by 



Hothorn et al. 



Herberich et al. 



(120081 ) and developed for the R software 



(120101 ). The max-t test does the simultaneous pairwise comparisons 
of means under control of the family-wise error rate, which is the probability of falsely 
rejecting the initial hypothesis (i.e., finding a significant difference among the means of any 
two groups in the data set even though there is actually no difference present). The test 
takes into account possible heteroscedasticity and the unequal sizes of the groups. 

We present the results of the tests under the form of simultaneous confidence intervals 
for all pairwise comparisons of group means. Confidence intervals including zero indicate 
no statistically significant differences. Confidence intervals near zero suggests some level of 
similarity. The farther from zero the more significant the differences. Note the smallness of 
some of the confidence intervals. This is due to the very large number of data used in each 
bin, which makes the tests extremely significant. 



In Figure [221 we show the confidence intervals associated with Figured (Section 13.11) . 
The SFG 2 galaxies, a) have more early-type, b) have older stellar populations, c) have 
more massive bulges, and d) are more luminous in B than the SFG 1. 



This manuscript was prepared with the A AS L^TjrjX macros v5.2. 
4 http://CRAN.R-project.org 
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In Figure [23] we present the confidence intervals associated with Figure [TT1 (Section 13. 2p 
where we compare in the SFG 1 the variations of physical parameters in 8 subsamples of 
metallicities: in a) the morphologies are significantly different except in pairs of bins (1,2) 
and (6,3). The general trend is for an increase of metallicity in early-type galaxies; in b) all 
the subsamples have different stellar population ages, the general trend being an increase of 
age with metallicity; in c) the variation in bulge mass is less obvious (many bins show no 
difference) but the trend is clear the mass increasing with the metallicity; in d) the increase 
in total mass with metallicity is more obvious. 

In Figure [2H we present the confidence intervals associated with Figure [T21 (Section 13. 2p 
where we compare in the SFG 2 the variations of physical parameters in 8 subsamples of 
metallicities. The same trends are observed than in the SFG 1. However, the variations 
are more subtile, which suggests that the parameters are somewhat independent of the 
metallicity; that is the SFG 2 are mostly early type galaxies (a) with massive bulges (c) 
dominated by old stellar populations, with a mass (d) increasing with the metallicity. 

In Figure we present the confidence intervals associated with Figure [H] (Section H|) 
where we compare the intensity of the bursts in the SFG 1 and SFG 2 in the two first bins 
in metallicities (1 and 5). The SFG 2 in bin 5 are clearly experiencing more intense bursts 
than the galaxies in all the other bins. However, also the SFG 1 in bin 5 (which are richer 
in nitrogen) are experiencing more intense bursts than the galaxies in bin 1 (poorer in 
nitrogen). In fact, the intensity of star formation in this subsample is comparable to what is 
observed in the SFG 2 in bin 1 (where the nitrogen abundance of the SFG 2 is comparable 
to that in the SFG 1) 

In Figure [26] we present the confidence intervals associated with Figure [17] (Section [4]) 
where we compare in the SFG 1 the variations of physical parameters in 9 subsamples 
increasing in redshift. The general trend observed in a) and b) is an increase of metallicity 
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and nitrogen at high redshift. There are fewer variations in morphological type (c) and 
stellar population (d), the trends being toward early-type but younger ages at high redshift. 
There is however a significant variation of bulge mass (e) and total mass (f), both increasing 
with the redshift. 

In Figure [27] we present the confidence intervals associated with Figure [TS] (Section H]) 
where we compare in the SFG 2 the variations of physical parameters in 9 subsamples 
increasing in redshift. Contrary to the SFG 1 we find almost no variation of metallicity 
in a) and nitrogen abundance in b) with the redshift. There are very few variations in 
morphological type (c) and stellar populations (d), the trends being toward late- type (the 
contrary than in the SFG 1) and younger ages at high redshift. However, like in the SFG 1 
there is a significant variation of bulge mass (e) and total mass (f), both increasing with the 
redshift. 
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Fig. 1. — Diagnostic diagram for the 122751 SDSS galaxies class ified as SFG by applyin g 
the empirical separation criterion (solid grey curve) proposed by iKauffmann et al.l ( 120031 ). 
According to this criterion, none of these galaxies should show an AGN contribution. 
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ig. 2 . — Same as Figure [T] using different separation criteria: a) according to 



Coziol et al. 



( 1999I ) models of HII regions, where Seq. 1 is a standard HII region abundance model, and 



Seq. 2 a similar model including an excess in nitrogen abundance by 0.2 dex-the median of 
these models (bold dashed curve) separates nitrogen poor from nitrogen rich S FGs; b) accord 



i ng to an alternate criterion to separate SFGs from AGNs, as proposed by 



Stasinska et al. 



(}2006|). 
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ig. 3. — Excitation test, to check for the presence of AGNs, as proposed by 



Coziol et al. 



( 1999I ). The lower left boxes contain galaxies with normal level of excitation (due to massive 



stars), while galaxies showing an excess of excitation in both lines fall out of the limits of 
these boxes. The four subsamples correspond t o nitr ogen poor SFGs (SFG 1) vs. nitrogen 



rich SFGs (SFG 2), according to 



Coziol et al. 



fjl999f ). and pure SFGs (SFG 3) vs. AGNs 



(SFG 4) according to Stasinska et al. (2006). 
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Fig. 4. — Second test for the presence of AGN based on the relations between the logarithm of 
the luminosity of the continuum (Lpc) and the logarithm of the luminosity of the ionized gas 
(L Ha ): a) results for the SFG 1 and SFG 2 samples; b) results for SFG 3 and SFG 4 samples. 
The significantly steeper lin ear relations for the Seyfert 2 and LINERs were determined by 



Torres-Papaqui et al. 



(120 121 ) . For numerical values see Table [U 
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Fig. 5. — Nitrogen abundance as a function of metallicity: in a) the SFGs poor in nitrogen 



(SFG 1), and in b) the SFGs r ich in nitrogen (SFG 2). A 



evolution models proposed by 



so shown are the closed-box chemical 



Vila- Cost as fc Edmunds! (119931 ): the secondary model (dashed 



line) and the primary + secondary model ( continuous curve). The two other curves trace 



the sequences obtained in 



Coziol et al. 



(119991 ) . In c) and d) we show the same results relative 



to the secondary relation for nitrogen enrichment, with their respective medians. 
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Fig. 6. — Mass-metallicity and mass-nitrogen excess relations for the SFG 1 and SFG 2. 
Overplotted on the data in a) an d b) we compare our linear correlations with the relation 
obtained by lTremonti et al.l ( 120041 ) . Results for the correlation tests are presented in Table EJ 
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Fig. 7. — Histograms for the metallicities and absolute magnitudes in the SFG 1 and SFG 2. 
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Fig. 8. — Box-whisker plots comparing the characteristics of the galaxies in the SFG 1 and 
SFG 2 samples: a) morphologies, b) mean stellar population ages, c) bulge masses, and c) 
absolute magnitudes in B. The points correspond to the means. Notches (barely visible due 
to the large size of the samples) that intersect suggest no significant difference in median. 
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Fig. 9. — Variation of abundances, where we have separated the whole sample in 4 sub- 
samples with increasing metallicity. In each subsample we further distinguish between low 
(subsamples 1, 2, 3, and 4) and high (subsamples 5, 6, 7, and 8) nitrogen abundances. 
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Fig. 10. — Variation of the star formation history (SFH) in the different subsamples defined 
in Figure [9j 
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Fig. 11. — Box- whisker plots for the 8 SFG 1 subsamples, as defined in Figure [9^l. We 
compare: a) the morphologies, T, b) the mean stellar ages, t stan c) the bulge masses, d) the 
absolute magnitudes in B. The diamond points with error bars correspond to the means. 
Notches that intersect suggest no significant difference in median. 
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Fig. 13. — Results of Starburst 99 models assuming instantaneous bursts with different 
IMF powers and upper mass limits. The data correspond to the medians and percentiles 
as measured in the subsamples in metallicity as defined in Figure [9] (identified by their 
numbers) . 



-48- 




Fig. 14. — Results of Starburst 99 models assuming continuous star formation with different 
IMF powers and upper mass limits. The data correspond to the medians and percentiles 
as measured in the subsamples in metallicity as defined in Figure [9] (identified by their 
numbers) . 
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Fig. 16. — Box-whisker plots for the burst intensities (in unit of M yrs 1 kpc 2 ) for the galax- 
ies in both samples showing the higher recent SFR (subsamples 1 a nd 5). The means with 



estima ted errors are also shown, as well as the threshold adopted by 



Strickland fc Heckman 



fl2009h . 
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Fig. 19. — SFH evolution with redshift in the SFG 1. The redshift of the galaxies increases 
from 1) to 9). 
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Fig. 20.— Same as Figured! for the SFG 2. 
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Fig. 21. — Illustration of the sequence of bursts model with starburst winds as main com- 
ponent. The star formation in the SFG 1 and SFG 2 can be reproduced by a sequence of 
bursts of star formation happening during their formation. Each burst causes the metallicity 
to increase and the ratio N/O to decrease. After 10 s yrs the intermediate mass stars formed 
during one burst begin to eject their nitrogen in the ISM, and the ratio N/O increases as the 
metallicity stay constant. The difference between the SFG 1 and SFG 2 is in the intensity 
of the bursts: in the SFG 2 more intense bursts, related to the formation of massive bulges, 
produce effective starburst winds, and the galaxies lose some of their oxygen. The number 
of intermediate stars being higher, the production of nitrogen is higher, and coupled with 
the lost in oxygen, a relative excess of nitrogen is observed for the same metallicity. 
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Fig. 22. — The confidence intervals associated with Figure El in Section I3TTT) : a) morphologies; 
b) stellar population ages; c) bulge masses; d) luminosity in B. 



-57- 



95% family-wise confidence level 



95% family-wise confidence level 



w 



-1.0 



-0.5 



o.o 



a) 

0l5 



b) 



-0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 



95% family-wise confidence level 



e confidence level 



HO 



£3 



( — •— i 



HH 

^02 ^Ol 00 Ol 02 olT 




-1.5 



-1.0 



^5 ocT 



d) 



0.5 



1.0 



Fig. 23. — The confidence intervals associated with FigurefTTlin Section l3~2l a) morphologies; 
b) stellar population ages; c) bulge masses; d) total mass. 
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Fig. 24. — The confidence intervals associated with FigurefT2~lin Section l3T2l a) morphologies; 
b) stellar population ages; c) bulge masses; d) total mass. 
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Fig. 25. — The confidence intervals associated with Figure [161 in Section HI 
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Fig. 26. — The confidence intervals associated with Figure IT71 in Section HJ a) metallicity b) 
nitrogen abundance; c) morphologies; d) stellar population ages; e) bulge masses; f) total 
mass. 
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Fig. 27. — The confidence intervals associated with Figure [18] in Section HJ a) metallicity b) 
nitrogen abundance; c) morphologies; d) stellar population ages; e) bulge masses; f) total 
mass. 
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Table 1: Slopes and correlation parameters for log L(Ha) vs. log L F c A4800. 



Sample 


slope 


Spearman 


Pearson 


SFG 1 


1.10 ±0.05 


0.883 


0.891 


SFG 2 


1.14 ±0.03 


0.932 


0.932 


SFG 3 


1.11 ±0.05 


0.877 


0.886 


SFG 4 


1.13 ±0.03 


0.931 


0.931 


Seyfert 2 


1.48 ±0.03 


0.810 


0.829 


LINER 


1.45 ±0.03 


0.739 


0.773 



Notes: the probability of correlation by chance is practically for all the tests. 
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Table 2: Correspondence between Hubble morphology and morphological indices, T 



Hubble 


E 


E/SO 


SO 


SO/Sa 


Sa 


Sab 


Sb 


Sbc 


Sc 


Scd 


Sd 


Sdm 


Im 


de Vaucoulcurs 


-6 to -4 


-3 to -1 







1.0 




3.0 




5.0 




7.0 


8.0 to 9.0 


10 to 11 


This work 


-5 


-2 





0.5 


1.0 


2.0 


3.0 


4.0 


5.0 


6.0 


7.0 


9.0 


10.0 
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Table 3: Pearson (r) and Spearman (r) correlation tests with chance probabilities (P) for 

the mass-metallicity and mass-nitrogen excess relations. 

Pairs r P(r) r s P(r s ) 

M B vs [0/H] (SFG 1) -0.4091 0.0450 -0.4369 0.0445 

M B vs [O/H] (SFG 2) -0.1376 0.4893 -0.1390 0.5143 

M B vs A(N/0) (SFG 1) -0.2016 0.3652 -0.2160 0.3865 

M B vs A(N/0) (SFG 2) -0.0638 0.5964 -0.0842 0.5643 
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Table 4: Physical characteristics of SFG 1 as divided by metallicity intervals 



Bin 


Number 


Range in 


Range in 


T 






M B 


id. 


of gal. 


log(0/H) 


log(N/0) 










1 


1759 


-3.85 < x < -3.55 


-1.85 < y < -1.40 


6.0 


7.88 


9.42 


-17.81 


2 


11375 


-3.55 < x < -3.25 


-1.60 < y < -1.15 


6.0 


8.25 


9.44 


-18.14 


3 


25561 


-3.25 < x < -2.95 


-1.25 < y < -0.80 


6.0 


8.44 


9.55 


-18.84 


4 


13697 


-2.95 < x < -2.65 


-1.05 < y < -0.60 


5.0 


8.53 


9.66 


-19.36 


5 


861 


-3.85 < x < -3.55 


-1.40 < y < -0.95 


6.0 


7.99 


9.46 


-18.24 


6 


5640 


-3.55 < x < -3.25 


-1.15 < y < -0.75 


6.0 


8.32 


9.51 


-18.66 


7 


5739 


-3.25 < x < -2.95 


-0.80 < y < -0.35 


4.0 


8.62 


9.75 


-19.42 


8 


3126 


-2.95 < x < -2.65 


-0.60 < y < -0.15 


4.0 


8.68 


9.77 


-19.47 
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Table 5: Physical characteristics of SFG 2 as divided by metallicity intervals 



Bin 


Number 


Range in 




Range in 




T 






M B 


id. 


of gal. 


log(0/H) 




log(N/0) 












1 


41 


-3.85 < x < - 


-3.55 


-1.85 < y < - 


-1.40 


6.0 


7.82 


9.52 


-18.23 


2 


207 


-3.55 < x < - 


■3.25 


-1.60 < y < - 


1.15 


5.0 


8.68 


9.57 


-18.08 


3 


9214 


-3.25 < x < - 


-2.95 


-1.25 < y < - 


0.80 


5.0 


8.53 


9.59 


-19.17 


4 


5544 


-2.95 < x < - 


-2.65 


-1.05 < y < - 


-0.60 


4.0 


8.62 


9.69 


-19.46 


5 


590 


-3.85 < x < - 


-3.55 


-1.40 < y < - 


0.95 


5.0 


7.87 


9.61 


-18.83 


6 


3327 


-3.55 < x < - 


-3.25 


-1.15 < y < - 


-0.75 


5.0 


8.47 


9.63 


-18.86 


7 


30201 


-3.25 < x < - 


-2.95 


-0.80 < y < - 


-0.35 


4.0 


8.71 


9.75 


-19.36 


8 


5863 


-2.95 < x < - 


-2.65 


-0.60 < y < - 


0.15 


4.0 


8.72 


9.78 


-19.42 
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Table 6: Intensity of star formation and absolute magnitude of SFGs showing evidence of 
recent bursts 



Sample SFR/Area M s 





Min 


1Q 


Median 


Mean 


3Q 


Max 


Min 


1Q 


Median 


Mean 


3Q 


Max 


SFG 1 (1) 


0.002 


0.020 


0.032 


0.040 


0.051 


0.351 


-22.79 


-18.39 


-17.72 


-17.81 


-17.19 


-16.02 


SFG 1 (5) 


0.002 


0.020 


0.036 


0.048 


0.062 


0.395 


-21.51 


-18.96 


-18.21 


-18.24 


-17.50 


-16.05 


SFG 2 (1) 


0.009 


0.024 


0.040 


0.054 


0.059 


0.257 


-20.42 


-18.81 


-18.23 


-18.23 


-17.64 


-16.44 


SFG 2 (5) 


0.002 


0.029 


0.062 


0.080 


0.105 


0.509 


-21.92 


-19.77 


-18.85 


-18.83 


-17.98 


-16.12 
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Table 7: Physical characteristics of SFG 1 as divided by redshift intervals 



Bin 


Range in z 


Number 


T 


L°g(^>L 


Log(^%f ) 


M B 


id. 




of gal. 










1 


0.025 < z < 0.050 


18778 


5.0 


8.40 


9.22 


-17.86 


2 


0.050 < z < 0.075 


19718 


5.0 


8.49 


9.47 


-18.68 


3 


0.075 < z < 0.100 


13479 


5.0 


8.43 


9.64 


-19.21 


4 


0.100 < z < 0.125 


8037 


5.0 


8.31 


9.82 


-19.67 


5 


0.125 < z < 0.150 


4682 


5.0 


8.26 


9.94 


-20.00 


6 


0.150 < z < 0.175 


1741 


5.0 


8.21 


10.08 


-20.34 


7 


0.175 < z < 0.200 


819 


5.0 


8.10 


10.20 


-20.63 


8 


0.200 < z < 0.225 


406 


5.0 


8.08 


10.26 


-20.82 


9 


0.225 < z < 0.250 


98 


5.0 


7.99 


10.31 


-20.93 
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Table 8: Physical characteristics of SFG 2 as divided by redshift intervals 



Bin 


Range in z 


Number 


T 


L°g(^>L 


Log(^%f ) 


M B 


id. 




of gal. 










1 


0.025 < z < 0.050 


7338 


4.0 


8.77 


9.35 


-18.26 


2 


0.050 < z < 0.075 


13063 


4.0 


8.83 


9.60 


-18.78 


3 


0.075 < z < 0.100 


12653 


4.0 


8.74 


9.78 


-19.21 


4 


0.100 < z < 0.125 


9609 


4.0 


8.58 


9.96 


-19.62 


5 


0.125 < z < 0.150 


6599 


4.0 


8.52 


10.08 


-19.98 


6 


0.150 < z < 0.175 


3008 


5.0 


8.37 


10.20 


-20.34 


7 


0.175 < z < 0.200 


1466 


5.0 


8.27 


10.33 


-20.63 


8 


0.200 < z < 0.225 


953 


5.0 


8.20 


10.40 


-20.83 


9 


0.225 < z < 0.250 


304 


5.0 


8.06 


10.47 


-21.02 



